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Abstract 



We present a systematic QCD analysis of the strange-charm and bottom-top con- 
tributions to transverse and longitudinal structure functions in charged-current deep 
inelastic scattering. Various 0{al) schemes are studied and compared. The depen- 
dence of their predictions on the factorization scale is investigated. The theoretical 
uncertainties resulting from the choice of the scheme and of /i^ are estimated. 



Neutrino Deep Inelastic Scattering (z/DIS) is an invaluable source of information on 
nucleon parton densities. In particular, neutrino reactions play a fundamental role in 
the determination of the strange quark distribution [|l|. 

Both to theorists and to experimentahsts neutrino-induced deep inelastic processes 
pose a number of subtle problems, due to the peculiar structure of the weak currents. 
These currents are not conserved, hence large longitudinal contributions are expected 
to exist , which badly break the Callan-Gross relation. Moreover, in charged-current 
weak DIS, flavours of different masses are mixed and interesting threshold effects occur 
0. For instance, the excitation of strangeness (a light flavour) produces charm (a 
heavy flavour) via the lowest order transition W~^s — > c, and there is the simultaneous 
production of s and c in the higher-order VT-gluon fusion reaction W~^g — sc. A 
similar situation takes place in the bottom-top sector. 

Two important questions arise in charged-current DIS. They possess a great theo- 
retical significance and, at the same time, are constantly faced by experimentalists in 
their analyses. 

The first question is whether a leading order (LO) approach, based on the quark 
excitation processes W~^q' q, W^q — ^ q' {q' = s,b; q = c,t), and on the slow rescal- 
ing mechanism (see below), is sufficient to carry out a reliable data analysis. The slow 
rescaling (SR) procedure is appealing because it is relatively easy to perform, but it is 
known to give unreliable results, in some instance. Previous theoretical |], H, |^, ^ 
and experimental studies on the determination of the strange density from neu- 
trino DIS have indeed taught us the importance of quark-mass corrections and current 
non-conservation effects, which manifest themselves through the next-to-leading order 
(NLO) VT-gluon fusion diagrams. However, when heavy quarks are involved, there is 
no unique way to take higher order effects into account in a partonic language. Various 
approaches exist, which hopefully (but not necessarily) are equivalent to each other. 

Here comes the second problem announced above: the treatment of heavy flavours. 
It seems natural to make the notion of "heavy" and "light" depend on the physical scale 
of the process (this has been emphasized in [^). At — 10 — 30 GeV^, which is the 
typical average scale of present neutrino experiments, s is light whereas c is undoubt- 
edly heavy. On the other hand, at ^ 10^ GeV^, where the bottom contribution 
becomes visible, s, c and even b can be safely considered light quarks and only t is really 
heavy. How should one treat "heavy" flavours in the large-Q^ limit, where collinear 
singularities emerge ? And how should one deal with the coexistence of "light" and 

^We shall refer to the 0{a'^) quark excitation processes as "leading order", and to the 0{al) W- 
gluon fusion processes as "next-to-leading order" . This terminology is admittedly ambiguous and is 
used here only for practical purposes and to adhere to some experimental papers 0| ; it differs from 
that adopted by other authors, e.g. M. 
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"heavy" quarks in the same QCD processes, which is typical of charged-current DIS 
and represents a further source of ambiguity ? 

Various schemes for a NLO QCD analysis of heavy quark contributions to elec- 
troweak structure functions are available on the market, but a detailed comparison of 
their predictions and their stability against changes of the factorization scale is still 
lacking. In this paper we want to bridge this gap, bringing to completion a program 
initiated in |igg. We shall offer a panoramic and detailed view of the heavy flavour 



sectors of F2 and F^, with their dependence on x, Q , and in particular on the NLO 



scheme and the factorization scale. Results for F3 have been presented elsewhere 

Let us start from the QCD factorization formula for the qq' contribution to the 
charged-current nucleon structure functions F2 and Fl (g = c, t; q' = s,b), that can 
be formally written as (0 means convolution) 

F^iiQ') = E ® ^2!l(/^', Q') ■ (1) 

a 

Here the sum is made over all parton species, Fj'^l ^^^^ structure functions for W 
scattering on parton a, and yU^ is the factorization scale. 

At order q;° only quarks and antiquarks contribute to the sum in (|l]) and Fg'*^, are 
proportional to the cross sections for the excitations W~^q' —>■ q, W^q g', which are 
essentially delta functions of the slow rescaling variables [0. For cs one has, e.g. 

F-(a;) = 2e[s(0+c(x)], (2) 

where ^ = (1 +Tn'^/Q'^) x. The longitudinal structure function, although nonvanishing, 
is suppressed by a factor ml/Q^ {rric is the mass of the charmed quark). Notice that 
eq. is a parton model formula: the only dependence is of a kinematical nature. 

The slow rescaling (SR) procedure consists in using eq. (||) with a dependence 
d la Altarelli-Parisi of the parton distributions. For cs 

SR : (x, Q') = 2^ Q') + c(x, Q^)] . (3) 

We recall that in some experimental papers f^, |^ a "leading order analysis" of data is 
performed, which is in fact the slow rescaling analysis sketched above. 

At order the main contribution to the qq' component of the structure functions 
is given by the VT-gluon fusion (GF) term, which reads 

OF : F^i{x,Q') = ^ /' ^ygiy,f.')Fl, (^,qA , (4) 

2tt Jax y \y I 



preprint has recently appeared Q which partly overlaps the study carried out in Ref. ||l^ and 
here. However a different perspective is adopted in 
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where a = 1 + (m^ + m'^)/(5^ and g{y-,iJ^) is the gluon density at the factorization 
scale. The exphcit expressions for and F| are, in the hmit m' — |13|, 
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gg splitting function expressed 



In eqs. s = Q'^{1 — z) / z and 'Pgg is the usual g 
in terms of the rescaled variable z = z{l + m^/Q^)- 

When the quark q' is a light flavour, which means either m' = or ^ m'^, 
the logarithm L{m', m) = log [{s — rn^Y / sm'"^] develops a (collinear) singularity. Same 
considerations apply also to the logarithm L{m,m') = log(,s/m^), when ^ m^. 
There are various methods to subtract (or to deal with) such collinear divergences. 

Let us focus for the moment on the charm-strange sector. We shall discuss three 
different prescriptions. Two of them are variations of the so-called "fixed flavour 
scheme" (FFS) [^]. In this scheme the charmed quark is not considered to be a parton, 
i.e. an internal constituent of the nucleon, and its 0{a^) excitation diagram is omitted. 
We call FFS (a) and FFS(b) the two realizations of the fixed flavour scheme that we 
shall consider in the following. They differ for the treatment of the strange quark: in 
FFS(a) strangeness is considered as a heavy flavour, with a mass ~ 0.3 — 0.5 GeV^; 
in FFS(b) the strange quark is viewed as a massless parton. The third prescription 
is the "variable flavour scheme" (VFS) in which charm is considered as a partonic 
constituent of the nucleon at large and /i^ ^ m^) and the 0{a^) charm 

excitation diagram is taken into account, with a proper subtraction term. 

In the FFS (a) scheme one treats the strange quark (but not the charmed one) 
as a parton, setting m' to zero and subtracting in eqs. (|^J^) the terms proportional 
to logQ^/m'^, which are embodied in the Altarelli-Parisi evolution of the strange 
distribution. In explicit form, the FFS (a) prescription reads 



FFS(a): F^Ux,Q') = 2^ s{^,Q') + 



27r 



dy 

ax y 



.y 



(7) 



■^For the sake of simplicity, we shall assume in our formulas that both the physical scale and 
the squared mass of the quark q are much larger than the squared mass m'^ of the quark q'; this 
is of course always true for the cs sector, and is true also for the bt sector at lO^GeV^. In the 

calculations the appropriate masses are used and no approximation is made. 
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where means eqs. (^H), with the terms proportional to logQ^/m'^ subtracted 
off. 

In the FFS(b) scheme eq. (^) is the only contribution to -^2'^^. A small but finite 
mass for the strange quark is retained and the full dependence on m and m' in the 
gluonic structure functions and (H). Of course no singularity appears. Explicitly 
one has 

FFS(b) : (x,Q^) = ^ f ^ y 9{y, f^') (-,qA , (8) 

Iir Jax y \y J 

with m = rric 7^ and m' = 7^ 0. 

In the fixed flavour schemes the number ^/ of active flavours is set to 3, irrespective 
of the scale of ct^. 

In the VFS scheme the strange quark is a massless parton and the charmed quark 
too is treated as a parton at large Q^. Both the strange and the charm 0{a^) excitation 
diagrams are considered but the dependence of the parton distributions is kept in 
the kernel of the 0{al) factorization integral. The collinear singularities are subtracted 
out by setting 

L(m,m') ^ log — , L(m',m)^ log , (9) 

where the subtraction scale is customarily taken to be equal to the factorization scale. 
Thus the VFS prescription is 

VFS: F£l(a:,Q^) = 2i[s{i, i?) + c{x, ix")] 

where means the gluonic structure functions, eqs. d^J^), with the two subtractions 
(P). Notice also that in this scheme the number of active flavours is 4. 

The bottom-top sector is treated along similar lines, but, out of the three ap- 
proaches that we have just outlined, only FFS(b) and VFS survive. 

What experimentalists do in their analyses is to use one of the methods presented 
above (including the slow rescaling prescription), choose a factorization scale and ex- 
tract the charm and strange (or bottom and top) densities from the measured structure 
functions. It is clear that there at least two sources of uncertainty in this procedure: 
one coming from the choice of the scheme, the other coming from the choice of /x^. 
The strategy adopted here reverses the experimental procedure: we use a set of parton 
distributions provided by a global fit and study the dependence on , /i^ and on the 
scheme of the charged-current structure functions in the heavy quark sector. 
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In our calculations we use the MRS- A fit of parton densities ||T5 . For consistency we 
adopt the same set of quark masses of Ref. that is = 2.7 GeV^, ml = 30 GeV^. 
We use also rrit = 180 GeV and, for the FFS(b), = 0.3 GeV. As for the factorization 
scale, authors who adhere to different NLO schemes give different suggestions: in the 
fixed flavour scheme of Ref. [|1^ , a value /i^ = is preferred, whereas in the variable 
flavour scheme of Ref. ||^ a value which asymptotically becomes of order of is 
suggested, and /i^ = is practically used. For homogeneity we decided to choose 

the same value for all schemes, thus we take /x^ = for the calculations presented 
in Figs. 1-4 and postpone the illustration of the /i^ dependence to Fig. 5. 

In Fig. 1 we present at two different values (Q^ = 25, 100 GeV^). The results of 
all prescriptions [FFS(a), FFS(b), VFS, SR] are shown. Consider, for instance, ^2'^'* at 

= 25 GeV^. If we take the difference between the results of the various prescriptions 
as a measure of the theoretical uncertainty arising from the choice of the scheme, this 
uncertainty amounts to ±3% at x = 10^'^, and to ±20% at x = 10^^. We have checked 
that the three NLO schemes give largely different predictions at lower x {x < 10~^); 
however a correct analysis of this region requires the use of the kt factorization method, 
which is beyond the purposes of this work. A forthcoming paper is specifically devoted 
to this issue |jl6|. Note that the slow rescaling curve almost coincides with the VFS 
result. On the other hand, the largest discrepancy occurs between the FFS(a) and the 
FFS(b) predictions, due to the 0{a^) term which is absent in FFS(b). We mention that 
the FFS(b) result for cs is affected by an additional (and large) uncertainty coming 
from the choice of m^. This makes FFS(b) intrinsically unreliable. 

If we confine ourselves to VFS and FFS(a), the theoretical uncertainty is ±(3 — 7)% 
over the whole x range at = 25 GeV^ and gets smaller with increasing Q^, being 
±(2 — 4)% at Q"^ = 100 GeV^. Thus it will be important to have data on F2'' at larger 

In Fig. 2 the curves for the longitudinal charm-strange structure function are 
displayed. The FFS(a) and FFS(b) results are much closer to each other in this case 
than in the F^^ case. The theoretical uncertainty amounts now to ±10% at = 25 
GeV^ over the whole x range. Obviously the SR curve for is much lower, being 
suppressed by a factor ml/Q^. The slow rescaling approach is expected to be unreliable 
when the longitudinal structure function dominates. The relevance of the longitudinal 
contribution is quantified by the ratio p = Fl/F2, which is a measure of the violation 
of the Callan-Gross relation. In the cs sector we find, in the FFS(a) scheme, p'^^ ~ 0.30 
at X = 10-3 ^ 10-2. 

The breaking of the Callan-Gross relation is stronger in the bt sector. In Fig. 3 
we show the SR, FFS(b) and VFS predictions for at two values {Q^ = 10^, 10^ 
GeV^). The difference between FFS and VFS amounts to ~ 10%, whereas the discrep- 
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ancy between slow rescaling and NLO schemes is more dramatic (up to ~ 70%). We 
do not plot F^*, limiting ourselves to mention that it dominates We find in fact 
p^\x = 10-3) = 0.97 {i.e. Fl* ~ F|*) at = 10^ GeV^, and p^*(x = 10-^) = 0.79 at 
= 10^ GeV^ (these values are almost scheme independent). 

For completeness we show in Fig. 4 the dependence of the structure functions 
^2^* and Ff'* in two schemes, VFS and FFS(a). The trend is similar in the two cases. 

Now, is there a way to prefer a scheme to another ? A practical criterion is to look 
at the dependence of the results on the factorization scale /i^. In Fig. 5 we present 
the /i^ dependence of at = 25 GeV^ in each scheme. It can be easily seen that 
VFS and FFS(a) are more stable against changes in /i^. A slight preference should be 
attributed to the latter scheme, which gives a practically constant prediction in the 
range /x^ = (5 — 100) GeV^, whereas in the same range the VFS result changes by 
~ 25% at X = 10^3, and by ~ 15% at x = 0.1. We have checked that also at = 100 
GeV^ FFS(a) is more stable than VFS in the range /i^ = (5 — 1000) GeV^. The main 
variation of VFS occurs for small /x^ values, so that if VFS is used a /x^ > is 
preferable. 

We remark here that in the next-to-leading CCFR analysis of dimuon data [^, 
which resorts to the VFS scheme, the factorization scale is chosen to be equal to twice 
the maximum available transverse momentum of the final charm quark, n = 2p™"^. At 
X = 10~^ this scale is about two orders of magnitude larger than the CCFR average 
Q^, which is 22 GeV^. The quoted uncertainty on the strange distribution due to the 
arbitrariness of the choice of /x^ is about 5% and is obtained by varying /x between 
ptnax r^pinax^ rcsults show that the CCFR Collaboration underestimates the 
theoretical uncertainties of its determination. 

Let us summarize our work. In this paper we provided an (x, Q^,/x^) map of the 
heavy quark contributions to charged current structure functions in some customary 
QCD schemes. The best way to choose the optimal scheme and factorization scale is 
to study the perturbative stability, which is possible only by pushing at least to order 

the knowledge of the factorization kernels. Some partial results heve been obtained 
in the electromagnetic case |l^, |l3, but a similar treatment of the weak DIS is still 
lacking and will be object of a future investigation. However some conclusions can be 
already drawn. While the slow rescaling mechanism should be definitely abandoned 
- at least when longitudinal contributions are present (it may nevertheless be useful 
for the analysis of purely transverse structure functions such as xF^ [0) -, it seems 
that, in the charm-strange sector, the fixed fiavour scheme with subtraction for the 
strange quark - the scheme we have called FFS(a) - is preferable due to its greater 
stability. In any case, the uncertainties coming from the choice of the scheme and of 
the factorization scale add up to some tens per cent, and require a better consideration 
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both on the experimental and on the phenomenological side. More data, especially at 
larger Q^, are needed in order to achieve a well-established knowledge of the strange 
and charm distributions. 
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Figure Captions 



Fig. 1 at = 25, 100 GeV^. The solid line corresponds to the VFS resuh [ see the 
text], the dot-dashed one to FFS(a), the dashed one to FFS(b), the dotted one 
to SR. 

Fig.2 Ff at = 25, 100 GeV^. The curves are labelled as in Fig. 1. 

Fig.3 F|* at Q2 = iqs^ GeV^. The curves are labelled as in Fig. 1. 

Fig.4 The dependence of Fl'* and Only the VFS prediction (solid curve) and 
the FFS(a) (dot-dashed) prediction are shown. 

Fig. 5 The dependence of on the factorization scale /j,"^ a,t x — 0.001, 0.01 and 0.1 
and Q'^ — 25 GeV^. The curves are labelled as in Fig. 1. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9601385vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9601385vl 
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